Updated data of the 153 Galactic globular clusters are 
used to readdress fundamental parameters of the Milky 
Way, such as the distance of the Sun to the Galactic center, 
bulge and halo structural parameters and cluster destruc- 
tion rates. We build a reduced sample, decontaminated of 
the clusters younger than 10 Gyr, those with retrograde 
orbits and/or evidence of relation to dwarf galaxies. The 
reduced sample contains 116 globular clusters that are 
tested whether formed in the primordial collapse. The 33 
metal-rich globular clusters ([Fe/H] > —0.75) of the re- 
duced sample extend basically to the Solar circle and dis- 
tribute over a region with projected axial-ratios typical 
of an oblate spheroidal, Ax : Ay : Az ps 1.0 : 0.9 : 0.4. 
Those outside this region appear to be related to accre- 
tion. The 81 metal-poor globular clusters span a nearly 
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from the central parts to the outer halo, although sev- 
q ■ eral clusters in the external region still require detailed 
studies to unravel their origin as accretion or collapse. A 
new estimate of the Sun's distance to the Galactic cen- 
ter based on symmetries of the spatial distribution of 116 
globular clusters is provided with an uncertainty consider- 
| ably smaller than in previous determinations using globu- 

■ lar clusters, Ro = 7.2±0.3kpc. The metal-rich and metal- 
OO | poor radial-density distributions flatten for Rqc < 2kpc 

■ and are well represented over the full Galactocentric dis- 
^ | \ tance range both by a power-law with a core-like term and 

I ■ Sersic's law; at large distances they fall off as ~ R^ 39 . 
iy) \ Both metallicity components appear to have a common 
1 origin, which is different from that of the dark matter 
"*^h \ halo. Structural similarities of the metal-rich and metal- 
O ^ poor radial distributions with the stellar halo are consis- 
' | tent with a scenario where part of the reduced sample was 
; I ■ formed in the primordial collapse, and part was accreted 
\ in an early period of merging. This applies to the bulge 

■ as well, suggesting an early merger affecting the central 
| parts of the Galaxy. The present decontamination proce- 
dure is not sensitive to all accretions (especially prograde) 
during the first Gyrs, since the observed radial density 

■ profiles still preserve traces of the earliest merger (s). We 
estimate that the present globular cluster population cor- 
responds to < 23±6% of the original one. The fact that the 
volume-density radial distributions of the metal-rich and 
metal-poor globular clusters of the reduced sample follow 
both a core-like power-law and Sersic's law indicates that 
we are dealing with spheroidal subsystems in all scales. 
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Abstract. 



1. Introduction 

Globular clusters (GCs) are potential witnesses of the for- 
mation processes that gave rise to the Milky Way. Because 
of their long-lived nature, GCs formed in the initial phases 
of the Galaxy may preserve in their structure and spatial 
distribution information that is essential to probe these 
early conditions. In this sense, derivation of the present- 
day spatial distribution of GCs as well as their physical 
properties can be used to infer on the Galaxy formation 
and evolution processes and better trace out the geometry 
of the Galaxy. 

Early models suggested that the Galaxy formed as a 
consequence of a monolithic, dissipative collapse of a single 
massive, nearly-spherical spinning gas cloud (e.g. Eggen, 
Lynden-Bell & Sandage 1962; Sandage 1990). Initial con- 
ditions of the collapse included low-metallicity gas and a 
nearly free-fall regime. This process should be reflected 
in the GC population as a homogeneity in certain pa- 
rameters, such as orbital motions and a restricted age 
range. However, later work presented observational evi- 
dence that the present-day GC population results not only 
from the primordial collapse but from merging and cap- 
tures of smaller neighbouring galaxies early in the his- 
tory of the Galaxy (Searle & Zinn 1978; Zinn 1980) or in 
more recent events such as the accretion and disruption 
of the Sagittarius dwarf spheroidal galaxy (Ibata, Irwin Sz 
Gilmorc 1994, Ibata et al. 11997(1 . On theoretical grounds, 
N-body simulations (using standard cosmological condi- 
tions such as cold dark matter) suggest that the hierar- 
chical merging of satellites might be the main building 
blocks of galaxy formation (Bcllazzini, Ferraro & Ibata 
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Present-day data picture the region interior to ~ 
lOkpc (inner halo/bulge) as formed essentially by the pri- 
mordial dissipative collapse (e.g. van den Bergh & Mackey 
2005), while the region external to ~ 15kpc (outer halo) 
was formed by later infall and capture of smaller frag- 
ments (Searle & Zinn lX3?£l Zinn lTMS) . Mackey & Gilmore 
PUUljl and Mackey & van den Bergh (|5UU5jl studied the 
properties of GCs by means of the horizontal-branch mor- 
phology, photometric and structural parameters. They 
found significant differences among age/metallicity sub- 
groups. Mackey & Gilmore (|2004[1 reported that ~ 30% 
of the Milky Way GCs have properties similar to those of 
the GCs in the LMC, SMC, Fornax and Sagittarius dwarf 
spheroidal galaxies. This suggests that a significant frac- 
tion of the Galactic GC population, in particular outer 
halo ones, has an extragalactic origin. 

Regardless of the origin, the Galaxy contains ~ 150 
GCs that are characterized by a bimodal metallicity dis- 
tribution, and distances to the Galactic center of up to 
~ 150 kpc. With respect to the metallicity vs. position re- 
lation, the metal-rich GCs used to be associated with a 
disk while the metal-poor ones with the halo (Kinman 
1X3551 Zinn 1X3551 Armandroff [X3J3). More recently the 
metal-rich GCs have been found to characterize a bulge 
population (Minniti 1X3351 Barbuy et al. 1X3331 Cote 1X3331 
van den Bergh [2000 ) . Chemical enrichment models of 
the Galaxy, especially for the central parts (Matteucci & 
Romano 1999.; Matteucci, Romano & Molaro 1999), pre- 
dict that the formation of the bulge occurred from the 
same gas, but even faster than the inner halo. 

Considerable efforts have been undertaken in the last 
~ 15 years to obtain fundamental parameters of globular 
clusters by means of accurate CCD colour-magnitude di- 
agrams e.g. for the central parts of the Galaxy (Barbuy, 
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Bica k Ortolani lT^gjl . Harris (1996), as updated 1 in 2003 
(and references therein), compiled parameters of GCs that 
we will adopt as the starting point for this work. Hereafter 
we refer to Harris' database as H03. 

Previous work have focused on the study of global 
properties and correlations among intrinsic parameters of 
the Galactic GC system, including as well the search for 
correlations with position in the Galaxy, e.g. Djorgovski 
& Meylan $Mty . van den Bergh (|2()03[> and Mackey & 
van den Bergh ijMTSjl . 

The significant improvement in the observational GC 
data is in itself a reason for checking if classically adopted 
values related to the early formation of the Galaxy and 
the GC system are still valid. In the present study we use 
an updated set of GC parameters, e.g. reddening, metal- 
licity, distance from the Sun, age and orbital motion to 
address their spatial distribution. Since one of the latest 
derivations of the Galactic center distance using GCs was 
made by Reid (1993 ) we also discuss that value now based 
on the updated GC data. 

Our basic approach is to decontaminate the GC sample 
of the clusters clearly not related to the primordial collapse 
of the Galaxy. We do this by identifying young GCs, those 
with retrograde orbits and/or related to the accretion of 
dwarf satellite galaxies. 

This paper is organized as follows. In Sect.|21we present 
the updated GC sample and isolate the clusters probably 
associated to the primordial collapse. In Sect.Olwe present 
projections of their positions onto the (x,y), (x,z) and (y,z) 
planes, derive the distance to the Galactic center and infer 
on Galactic structure. In Sect. 0] we build GC radial den- 
sity profiles, fit them with different analytical functions 
and discuss GC destruction rates. In Sect. we discuss a 
possible scenario to account for the present spatial distri- 
bution of GCs. Concluding remarks are given in Sect. El 

2. The updated globular cluster data set 

In the last 10 years new entries have been added to 
the census of Galactic GCs either by means of discov- 
eries or identifications of misclassified open clusters, e.g. 
IC1257 (Harri s et a l. RWi\ . ESO280-SC06 (Ortolani, 
Bica & Barbuv EUUUjl . 2MASS-GC01 and GC02 (Hurt et 
al. I2UUU|) . GLIMPSE-C01 (Kobulnicky et al. l2TjU5|> . the 
diffus e clus ter-type object SDSSJ1049 + 5103 (Willman 
et al. I2()05|) . and the recently re-classified GC Whiting 1 
(Carraro 2005). For new and already known GCs H03 has 
provided constant updating of fundamental parameters. 

Parameters of the 153 presently known GCs in the 
Galaxy are listed in Tabled We complemented the H03 
database by adding information on new GCs in the last 2 
years as indicated in the notes to Tabled The presently 
updated GC data set will be hereafter referred to as the 
GC05 sample. 

To estimate errors in distance determinations for the 
subsequent analyses we took into account GC05 and typ- 

1 at http://physun.physics.mcmaster.ca/Globular.html 



ical distance errors in the literature which are dominated 
by reddening uncertainties (e.g. Barbuy, Bica & Ortolani 
1998). We adopt the following values for the reddening 
and distance error relation: e = 0.05 for E(B — V) < 0.2, 
e = 0.10 for 0.2 < E(B-V) < 0.8, e = 0.15 for 
0.8 < E(B - V) < 1.5, e = 0.2 for 1.5 < E(B - V) < 2.0, 
and £ = 0.25 for E(B — V) > 2.0, where £ is the fractional 
error in distance. Uncertainties in related parameters are 
obtained by propagation of £. 

As discussed in Sect. ^ the Galactic GC system is 
not expected to be homogeneous in terms of cluster ori- 
gin. Since we intend to base our analysis on GCs with a 
high probability of formation in the primordial collapse, 
we exclude from GC05 those with evidence either of ex- 
tragalactic origin or formation later than the collapse. To 
this category belong the GCs with retrograde orbits (e.g. 
Dincscu, Girard & van Altena H999|l . ages younger than 
lOGyr (e.g. Salaris & Weiss 121021 ■ direct relation or tidal 
debris of dwarf galaxies (e.g. Forbes, Strader & Brodie 
2004 ) and finally, luminous GCs with evidence of being ac- 
creted dwarf galaxy nuclei (e.g. Mackey & van den Bergh 
2005 ). References are given in col. 12 of Tabled We found 
that 37 GCs (24% of GC05) fit into one or more of these 
categories. The remaining 116 GCs are part of what we 
define hereafter as the reduced sample (RS-GC05). One 
caveat with respect to the decontamination process is that 
it certainly is not sensitive to all events, accretions in par- 
ticular, dating back to the first Gyrs of the Galaxy. For 
instance, extragalactic GCs accreted in this period with 
prograde orbits would in the present hardly be distinguish- 
able from the Milky Way's native population. 

Tabled contains, by columns: (1) - main GC designa- 
tion; (2) and (3) - Galactic coordinates; (4) - reddening; 
(5) - metallicity; (6) - distance from the Sun; (7) - input 
Galactocentric distance; (8) - output Galactocentric dis- 
tance fSect. ETl")) : (9)-(ll) - input Heliocentric components 
(Sect. |3J); (12) - indicators of non-collapse membership; 
and (13) - additional GC designations as compiled from 
the literature by one of us (E.B.). 

In Fig. n (upper panels) we compare the globular clus- 
ters of RS-GC05 with those in GC05 in terms of metal- 
licity. The well-known bimodal distribution in metallic- 
ity (seen e.g. as early as in Zinn|l985) is confirmed not 
only in the present GC05 sample (panel (a)) but also 
in RS-GC05 (panel (b)), however with a smaller ampli- 
tude ratio between the metal-rich and metal-poor GCs. In 
the subsequent analysis we adopt [Fe/H] = —0.75 as the 
metallicity threshold between metal-rich and metal-poor 2 
GCs. RS-GC05 contains 81 metal-poor GCs, 33 metal- 
rich and 2 with unknown metallicity (2MASS-GC02 and 
GLIMPSE-C01). The reddening distribution of the RS- 
GC05 GCs (panel (c)) is compared to those of the cor- 
responding metal-poor (panel (d)) and metal-rich GCs 
(panel (e)). The reddening distribution of the metal-poor 
GCs presents a maximum around E(B — V) « 0.05 and 

2 For simplicity we refer as metal-poor GCs the genuine to- 
gether with the intermediate-metallicity ones. 
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Fig. 1. GC metallicity distribution of the GC05 (panel 
(a)) and reduced samples (panel (b)). The adopted thresh- 
old between metal-rich and metal-poor GCs is indicated. 
The reddening distribution of the reduced sample GCs 
is in panel (c), the corresponding metal-poor ones are in 
panel (d) and the metal-rich ones in panel (e) . The insets 
in panels (c) - (e) show the high-reddening range. Note 
that the two GCs with E(B - V) « 4.8 and 5.5 (inset of 
panel c) have no metallicity determinations. 



a smaller one at E(B — V) « 0.45. The low-reddening 
values are related mostly to halo GCs, while the high- 
reddening ones belong to a more central metal-poor com- 
ponent that spatially coexists with the metal-rich bulge 
clusters (Barbuy, Bica & Ortolani 1998). The metal-rich 
GCs, in contrast, present a rather uniform distribution in 
the range 0.05 < E(B - V) < 1.3. 

3. Galactic structure and the distance to the 
center 

Inferences on the geometry of the GC system can be made 
by means of cluster positions projected onto the (x,y,z)- 
heliocentric coordinate axes (Table In this coordinate 
system the x-direction increases from the Sun towards the 
Galactic center, y is positive for £ = 0° — 180° and z in- 
creases towards the north Galactic pole. We consider sep- 
arately the metal-rich and metal-poor GCs, both of the 
GC05 and RS-GC05 samples. In Fig. [3 we show the po- 



sitions of the metal-rich GCs projected onto the (x,y), 
(x,z) and (y,z) planes for the GC05 (left panels) and RS- 
GC05 (right panels) samples. Because the (x,y,z) coordi- 
nates are in the heliocentric system, the centroids of the y 
and z distributions coincide with the Galactic center, while 
that of the x-coordinate is shifted from x = fSect. f5~T|) . 
GC05 metal-poor and metal-rich GCs are more widely dis- 
tributed than those in RS-GC05. This effect is minimized 
in the RS-GC05 plots because most of the GC05 outliers 
belong e.g. to accreted dwarfs or their debris, have young 
ages and/or retrograde orbits (Table The 2 remaining 
outliers in the metal-rich RS-GC05 (panels (b), (d) and 
(f)), NGC 6356 and Palomar 11, basically define the outer 
limits of the metal-rich system, slightly beyond the present 
determination of the Solar circle (see below) - 8 kpc (Reid 
1993}. These clusters deserve further attention to clarify 
whether they are young GCs, thus not related to the col- 
lapse, and/or located in the apogalacticon of their orbits. 

From the plots involving the RS-GC05 sample (Fig. [2J 
we estimate that the metal-rich GCs distribute essentially 
in a region with dimensions Ax ~ 12 kpc, Ay ~ 11 kpc 
and Az ~ 5 kpc, corresponding to axial-ratios Ax : Ay : 
Az k 1.0 : 0.9 : 0.4. These axial-ratios can be accounted 
for by an oblate spheroid with ~ 5.5 kpc in radius and 
~ 2.5 kpc in height, a structure spatially coincident with 
the bulge. As compared to earlier studies the present x 
and y distributions of the metal-rich GCs are similar in 
extent (Fig. |2J) because of the minimization of observa- 
tional errors achieved with present-day data. 

The metal-poor GCs of RS-GC05 are essentially con- 
tained in a region with dimensions Ax ~ 35 kpc, Ay ~ 
36 kpc and Az ~ 30 kpc, with axial-ratios s» 1.0 : 1.0 : 0.8. 
Considering uncertainties these axial ratios describe a 
slightly flattened sphere that reaches into the outer halo. 

We also infer on the spatial distribution of the 
GCs in RS-GC05 by means of the distribution func- 
dN 

tion 4>(0 = j which counts the number of GCs in 

bins of A£ = 1 kpc, for the x, y and z coordinates. 
According to the definition, 0(£) is related to the pro- 
jected one-dimensional number-density of GCs along a 
given direction. Fig. ^ shows the distribution functions 
of all GCs in RS-GC05, and the corresponding metal- 
poor and metal-rich ones, separately. As expected from 
Figs. [5] and |3 the distribution functions in y and z are 
symmetrical with respect to the centroid of the coordi- 
nate system (the Galactic center), while the shift in x 
provides the distance of the Sun to the Galactic center 
fSect. l3~T|l . The distribution functions in Fig.0|can be fit- 
ted both with exponential-decay and squared-hyperbolic 
secant functions. Exponential-decay functions usually de- 
scribe projected surface-density profiles in spiral galaxy 
disks (Binney & Tremaine I1987[) . while self-gravitating 
isothermal models such as the squared-hyperbolic secant 
have been used in edge-on disks (e.g. Rice et al. I1996f) 
and lenticular galaxies (van der Kruit & Searle in)81(l . Our 



purpose in fitting the distribution functions with a sym- 
metrical profile is to derive the distance of the Sun to the 
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Fig. 2. Spatial projections of the heliocentric positions 
of the metal-rich GCs of the GC05 (left panels) and re- 
duced (right panels) samples. The reduced sample pro- 
duces more concentrated distributions. The outlier metal- 
rich GCs NGC6356 and Palomarll are identified in pan- 
els (b), (d) and (f). 



Galactic center (Sect. ETT} . In this sense we adopted as fit 

I t-ta 

the exponential-decay </>(£) oc e I 4h function, since the 
respective correlation coefficients resulted larger than with 
squared- hyperbolic secant functions. The fits are shown in 
Fig-E and the resulting central positions (£ Q ) and scale- 
lengths are given in Tabled The distribution peaks in 
y and z occur, within uncertainties, at y = z = (Table 2 
and Fig.|2J|. 

Except for the central point in the y-distributions of 
the RS-GC05 sample (panel (b)) and corresponding metal- 



Fig. 3. Same as Fig.Elfor the metal-poor GCs. 



poor GCs (panel(e)) the exponential-decay function ac- 
ceptably fits the observed data for distances of up to 
~ ±7 kpc with respect to the peak, within uncertainties. 
At such distances we are probing not only the bulge but 
the inner halo as well. The fits preserve the symmetrical 
character of the observed profiles, and do not affect the 
centroid determination. This suggests that not many GCs 
remain undetected towards the central parts, at least to 
the point of affecting the distance determination. 

The individual fit of an exponential-decay profile to 
each of the (x,y,z) components does not necessarily imply 
a disk structure for the GC subsystems, since we are deal- 
ing with one-dimensional distribution functions and not 
surface density profiles. 

The scale-length ratios (Table I2J derived from the 
exponential-decay fits agree, within uncertainties, with the 
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Fig. 4. One dimensional distribution functions of the re- 
duced sample (top panels), metal-poor (middle panels) 
and metal-rich GCs (bottom panels). The profiles were fit- 
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Table 2. Parameters of the one-dimensional distributions 



Reduced Sample GCs 



(1) 


All 

(2) 


Metal-Poor 
(3) 


Metal-Rich 

(4) 


x (kpc) 


7.2 ±0.3 


7.1 ±0.5 


7.7 ±0.3 


x h (kpc) 


2.9 ±0.4 


2.9 ±0.6 


2.7 ±0.4 


y (kpc) 


0.0 ±0.3 


0.0 ±0.4 


0.1 ±0.3 


yh (kpc) 


1.9 ±0.4 


2.8 ±0.5 


2.1 ±0.4 


z (kpc) 


0.1 ±0.2 


0.4 ±0.3 


0.1 ±0.1 


Zh (kpc) 


1.8 ±0.2 


2.7 ±0.5 


1.1 ±0.1 


xh/yh 


1.5 ±0.4 


1.0 ±0.3 


1.3 ±0.3 


Xh/Zh 


1.6 ±0.3 


1.1 ±0.3 


2.4 ±0.4 


yh/zh 


1.1 ±0.3 


1.0 ± 0.2 


1.9 ±0.4 


Table Notes. Parameters of the function 0(£) 
to the £ = (x, y, z) distribution functions. 


oc e 1 £ h fitted 
Col. (2): All GCs 



of the reduced sample. 



axial-ratios estimated from Figs. [21 and |3 both for the 
metal-poor and metal-rich GCs of RS-GC05. 



3.1. Distance to the Galactic center 

The distance of the Sun to the Galactic center (Ro) has 
been a recurrent topic in the literature since Shapley's 
attempt in 1918 to derive it with globular clusters that 
resulted in Ro = 13 kpc. Since then different methods 
with more accurate data and larger GC samples have been 
used fot the same purpose. For instance Frenk & White 
(1982) using a sample of 65 metal-poor and 11 metal- 
rich GCs limited in latitude to avoid exceedingly large 
reddening errors affecting distances derived Ro = 6.8 ± 
0.8 kpc. 

Reid l|1993[l reviewed several estimators to derive Ro, 
among them the available GC parameters at that time. 
Estimates based on those GCs put Ro in the range 6.2- 
10.1 kpc. The average Ro from the GC determinations in 
Table 2 of Reid l|1993[) is 7. 9 ±1.4 kpc, which coincides with 
his best value of 8.0 ±0.5 kpc considering all methods, e.g. 
calibration by OB stars and Hi and Hn regions, GCs, RR 
Lyrae and red giants, among others. Since then this value 
has been widely employed in the literature. However, from 
his Fig. 3 it is clear that the x coordinates of the available 
GCs suffer from reddening/distance uncertainty effects. 
With the GC samples available at the time Maciel ([1993) 
and Rastorguev et al. I|1994|l obtained Ro ~ 7.6 kpc and 
Ro ~ 7.0 kpc, respectively. 

More recently, Eisenhauer et al. ( 2003) used VLT spec- 
troscopic observations of the orbit of the star S2 around 
SgrA* (assumed to be at the very center of the Galaxy) to 
derive Ro = 8.0 ± 0.4 kpc. In the same work they provide 
as well the value of Ro based on the statistical parallax 
distance of 106 late-type and 27 early-type stars located 
in the central 0.5 pc. They obtained Ro = 7.2 ± 0.9 kpc. 

One caveat is that the total-to-selective absorption ra- 
tio Ry = Ay/E(B — V) is not expected to be uniform. 
Variations of Ry in different directions throughout the 
Galaxy can occur (e.g. Sumi 120041 Ducati, Ribeiro & 
Rembold 2003J) • R-v is also affected by the effective wave- 
length shift in the filters owing to metallicity differences 
and reddening amount (Barbuy, Bica & Ortolani [1998 
and references therein). Detailed analyses of Ry in the 
directions of all GCs would be necessary to minimize Ry- 
related uncertainties. However, to a first approximation 
we assumed the H03 distances in Tabled H03 took into 
account the metallicity dependence of the absolute mag- 
nitude of the horizontal branch and from their data it can 
be inferred that a constant value of Ry = 3.1 was adopted 
throughout. 

At the 1(7 level the values of Ro provided by the 
one-dimensional exponential-decay fits of the metal-poor 
(panel (d) of Fig. 0J) and metal-rich GCs (panel (g)) are 
basically the same (Table 2). In this sense, to increase the 
statistical significance of the determination we applied the 
fit to the 116 GCs of RS-GC05 (panel (a) of Fig.©. We 
obtained an average value of Ro = 7.2 ± 0.3 kpc. This 
value puts the Sun « 0.8 kpc closer to the Galactic center 
than either the best one adopted by Reid i|1993|l or that 
derived by Eisenhauer et al. {2003). However, the present 
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value coincides with that of central stars by Eisenhauer et 
al. (HUH)). 

The present determination is based on a more accurate 
and numerous GC database, and consequently, the uncer- 
tainties in the value of Ro are a factor of ~ 3 smaller 
than in previous studies using GCs (see Table 2 of Reid 
1993). We used the new Rq determination to recalculate 
the Galactocentric distances (col. 8 of Table QJ. 



to avoid undersampling with increasing Galactocentric 
distance. Taken as face value the radial distribution func- 
tion as defined above should be applied to spherically 
symmetric systems, which is not the case of the oblate 
geometry of the metal-rich GCs of RS-GC05 (Sect. EJ). 
Implications of this difference in geometry will be dis- 
cussed in Sect. 14.41 



3.1.1. Variable total-to-selective absorption 



Following the analysis of Barbuy, Bica & Ortolani (1998 ) 
of the central Galaxy, we now explore the effect of a vary- 
ing total-to-selective absorption as a function of metal- 
licity and reddening amount for the whole GC system. 
According to Grebel & Roberts l|1995f) we adopted Ry = 

3.6 for the GCs more metal-rich than solar metallic- 
ity, and Ry = 3.1 for [Fe/H] < —1.0. Interpolation 
is used for intermediate values of metallicity. To the 
metallicity-interpolated Ry we add a further correction re- 
lated to reddening, AR V = 0.05 x E(B - V) ( Olson IHT75|) . 
Dependence of distance on varying Ry can be expressed 

(Ry-R^)E(B-V) 

as Ro = Ro x 10 b 

We applied the above corrections to the data in Tabled 
for all metal-rich GCs of RS-GC05 individually, leading to 
a smaller value of Ro = 6.6±0.5 kpc. Applying the same to 
all metal-poor GCs of RS-GC05 individually the distance 
of the Sun to the Galactic center remains essentially the 
same as before Ro = 7.3 ± 0.5 kpc (Table GJ. 

Finally, we consider the ensemble of the metal-rich 
GCs (Fig.^l in order to minimize individual uncertainties. 
The average metallicity and reddening are [Fe/H] as —0.55 
and E(B — V) ~ 0.5, providing an average R v = 3.35. 
For the metal-rich GCs Ry = 3.1 (Sect. EHJ and Ro = 

7.7 ± 0.5 kpc (Table |3J, and the resulting distance is 
Ro' = 7.3 ± 0.3 kpc, thus fully compatible with Ro de- 
rived from the metal-poor GCs. 

Irrespective of the metallicity and reddening law vari- 
ations for the metal- rich GCs, the present distance of the 
Sun to the Galactic center determination Ro « 7.2 kpc is 
robust, since it depends essentially on the larger sample 
of metal-poor GCs. This is due to the metal-poor GCs 
being rather insensitive to Ry assumptions and the fact 
that the current accuracy on their distances is significantly 
improved. 

4. Radial distribution of globular clusters 

The distribution in Galactocentric distance of the GC 
, , dN dN 

number-density, p(R) = - — = ttttj ls a potential 

dv 47rR^dR 
source of information not only on the present-day Galactic 

structure but the formation processes as well. To investi- 
gate this we build radial distribution profiles for the metal- 
rich and metal-poor GCs of RS-GC05 separately. Bins in 
radius of AR = 0.5 kpc for Rgc < 10 kpc are used to 
better sample the inner regions, while AR = 2 kpc for 
10 < Rcc(kpc) < 20 and AR = 10 kpc for R GC > 20 kpc 



4.1. Metal-poor GCs 



The radial distribution of the metal-poor GCs of RS-GC05 
is shown in panel (a) of Fig. A fraction of 74% of the 
81 metal-poor GCs is located at Galactocentric distances 
Rqc < 10 kpc, and 20% at R GC > 15 kpc (outer halo). 
The distribution falls off smoothly as a rather steep power- 
law ~ R-(3-6±o.2) for Rgc > 3.5 kp C . However, it flattens 
out for smaller Galactocentric distances as ~ R j -( 1 - 7 ± - 3 ). 
At least part of the flattening might be attributed to 
completeness effects in the crowded central region of the 
Galaxy. However, a near-IR survey with 2MASS (Dutra 
& Bica 2000) did not reveal any new GC in the central 
region. Two recent GC discoveries with 2MASS are not 
centrally located, since they are at I w 10° and near the 
plane, about halfway from the Sun to the Galactic center 
(Hurt et al. [2000 ) . Alternatively, the flattening for small 
Rgc may result from the cumulative destruction of GCs 
close to the Galactic center over a Hubble time (a pro- 
cess that certainly played a major role in depleting the 
original GC population - see Sect. 14. 5f) . or it may be an 
intrinsic feature of the radial distribution. Modelling of 
the spatial distribution of the old halo GCs beginning at 
the primordial collapse with cold baryonic gas and dark 
matter conditions suggests that the inner flattening may 
result not only from tidal destruction, but may in part be 
of primordial origin (Parmentier & Grebel [2005). 

Because of the flattening at small Rgc the simplest 
fits of the observed radial density profile are obtained 
with analytical functions that contain a core-like term. 
Following Djorgovski & Meylan (1903} we employ the 
function p(R) = Po/(l + R/Rc) Q , where Rc is the core- 
like radius. We will refer to this function as the composed 
power-law. The agreement between fit and observed radial 
distribution along the full Galactocentric distance range is 
excellent (Fig. 03 panel (a)). The resulting parameters are 
p = 3.9±2.5kpc- 3 , R c = 1.5±0.6kpc and a = 3.9±0.3, 
with a correlation coefficient CC = 0.88. 

Alternatively, in the inset of panel (a) we fit the 
metal-poor observed radial profile with Sersic's (|19rifi|l law, 
p(R) — ae _b [^ R °' 1 >_1 ] . Since it is rather insensitive to 
variations of Rc, we used the same core-like radius as that 
indicated by the composed power-law in order to have less 
free parameters when fitting Sersic's law. The best fit was 
obtained with n = 4.1 ± 0.7 and CC = 0.88. 
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4.2. Metal-rich GCs 

The 33 metal-rich GCs of RS-GC05 are contained in the 
region 0.66 < Rcc(kpc) < 8.3 (panel (b)), which shows 
that a sharp radial cutoff thus occurs in the metal-rich 
distribution near the Solar circle. Metal-rich GCs in GC05 
located outside this region appear to be related to accre- 
tion of dwarfs and/or young ages (Table ^). This con- 
trasts with the metal-poor GCs that distribute in the 
range 0.36 < Rgc (kpc) < 123. Similarly to the metal- 
poor GCs, a flattening in the radial distribution of the 
metal-rich GCs with respect to the extrapolation of the 
large Galactocentric distance power-law ~ R-( 3 - 2 ±°- 2 ) oc- 
curs for Rqc < 2kpc. This effect should be expected, 
since there is a lack of correlation of metallicity and GC 
luminosity (e.g. Djorgovski & Mevlan ll994l van den Bergh 
2003). Parameters of the composed power-law fit are p Q = 
6.3 ± 5.2kpc~ 3 , R c = 1.2 ± l.Okpc and a = 3.9 ± 1.2, 
with CC = 0.88. Sersic's law (inset of panel (b)) provides 
a fit with the exponent n = 2.9 ± 1.5 (CC = 0.87). In this 
fit we used Rc = 1.2 ± l.Okpc, as indicated by the com- 
posed power-law. The observed distribution (panel (b) of 
Fig. |3J cannot be fitted with an exponential-decay law, 
which precludes the presence of a disk. 

Probably as a consequence of the bulge/halo transi- 
tion, the flattening in both metal-poor and metal-rich ra- 
dial distributions begin at Galactocentric distances com- 
patible with the dimension of the bulge (Sect.|3J), partic- 
ularly with the (x,y,z) scale-lengths of the metal-rich GCs 
(Table EJ). 

Despite the marked difference in the radial extent of 
the metal-rich and metal-poor GC profiles, both distribu- 
tions present similar structural features such as flattening 
in the central region, core-like radius (Rc = 1.2 — 1.5 kpc), 
composed (a = 3.9) and single power-law slopes (n = 
3.2 — 3.5) and Sersic's law index (n = 2.9 — 4.1), within 
uncertainties. These similarities suggest that most of the 
GCs in both metallicity classes share a common origin. 



4.3. All GCs of the reduced sample 

The best-fit of the composed power-law to the radial dis- 
tribution of the 116 RS-GC05 GCs was obtained with 
p = 7.2±3.0kpc- 3 , R c = 1.9±0.5kpc and a = 4.4±0.3, 
with CC = 0.91 (panel (c) of Fig. EJ. Because the metal- 
rich GCs are contained in the region Rgc (kpc) < 8.3 
the slope of the composed power-law resulted slightly 
steeper than those of the metal-rich and metal-poor dis- 
tributions, as expected. In addition, the single power-law 
extrapolation for Rgc > 3.5 kpc falls off as ~ g_— (3.9±o.i)^ 
which within uncertainties basically agrees with that of 
the metal-poor GCs. The best Sersic's law fit was obtained 
with n = 3.0 ± 0.4 and CC = 0.89. Qualitatively, Sersic's 
law and the composed power-law provide essentially the 
same fit to the observed radial profile. 




1 10 100 

Rgc (^C) 



Fig. 5. Radial density profiles of the GCs in the re- 
duced sample as a function of Galactocentric distance. 
Panel (a) - metal-poor GCs; Panel (b) - metal-rich; 
Panel (c) - all GCs. Dashed line: single power-law fit for 
large Galactocentric distances. Solid line: fit of /?(R) = 
Po/(l + R/R C ) Q - Insets: fit of Sersic's law p(R) = 

ae -b[(R/R c )< 1 /"»-l]_ 

4.4. Spherically symmetric-volume densities for the 
oblate metal-rich sub-system 

For practical purposes we assumed spherical symmetry 
in the above analysis of metal-rich and metal-poor ra- 
dial density profiles. However, the metal-rich GCs dis- 
tribute through a region whose geometry is clearly oblate 
(Sect. |3J). Thus, spherically symmetric- volume densities 
calculated in radial bins beyond a few kpc from the center 
will be artificially decreased, as compared to those mea- 
sured in a genuinely (or approximately) spherical system. 
Consequently, both the measured power-law fall-off at 
large radii and central flattening degree in the metal-rich 
sub-system might result enhanced relative to the nearly 
spherically symmetric metal-poor sub-system. 

To investigate the effects of non-sphericity in the 
metal-rich sub-system we apply a coordinate transforma- 
tion to correct for its oblateness, y — > y/0.9 and z — » z/0.4 
(Sect.[3J). Subsequently we recalculate Galactocentric dis- 
tances, Rgc = -\/x 2 + y 2 + z 2 , and volume densities, 
p(R). Compared to the observed oblate profile (bottom 
panel) the corrected one (top panel of Fig. |SJ) presents 
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Fig. 6. Radial density profiles of the metal-rich GCs in 
the reduced sample measured with spherically symmetric 
volume densities on oblateness-corrected (top panel) and 
oblate (bottom panel) spatial geometries. 



similar shape and somewhat scaled-up Galactocentric dis- 
tances. Fit of the composed power-law results in a core- like 
R c = 0.9±0.9kpc and slope a = 3.2±0.9. Both R c and a 
are similar to the previous ones, but the slope is somewhat 
natter, as expected. 

Within uncertainties we conclude that measuring 
spherically symmetric volume densities in the oblate 
metal-rich GC sub-system has a small effect on structural 
parameters such as power-law fall-off at large radii and 
flattening degree in the central region. The effect is mini- 
mal probably because of the relatively small radial exten- 
sion (Rqc ~ 10 kpc) of the metal-rich sub-system. This 
effect is negligible for the nearly-spherical metal-poor GC 
sub-system (Sect.OJ). 

4.5. Globular cluster destruction rates 

The discussions in the previous section raised the question 
whether the flattening observed in the radial distribution 
for regions interior to Rgc ~ 2 kpc is a primordial feature 
or a consequence of enhanced GC-destruction rates near 
the Galactic center. Although the present analysis does 
not answer this question, it can be used to provide an 
estimate on the fraction of primordial GCs still present in 
the Galaxy. 



The Galactic environment, particularly near the cen- 
ter, tends to destroy star clusters because of enhanced 
tidal truncation and gravitational shocks due to passages 
close to the bulge and through the disk. The bulge, in 
particular, is very efficient in destroying clusters on highly 
elongated orbits (Gnedin & Ostriker ll997f) . and the pres- 
ence of the bar increases the destruction rates by providing 
a means to bring more clusters close to the Galactic center 
(Long, Ostriker & Aguilar 1992). Dynamical evaporation 
is probably the most important destruction mechanism in 
the present-day Galactic environment (Gnedin & Ostriker 
nW7j) . Aguilar, Hut & Ostriker l|1988|l estimate a current 
GC depletion rate of ~ 5% due basically to dynamical 
evaporation, indicating that most of the destruction took 
place in the past with bulge-shocking as the main factor. 
Hut & Djorgovski (1992) estimate that the present GC 
evaporation rate may be 5 ± 3 Gyr" 1 . 

Gnedin & Ostriker (1997) found larger destruction 
rates than previous work, predicting that 52% - 86% of 
the present GC population may be destroyed in the next 
Hubble time. They conclude that the present GC popula- 
tion must be a small fraction of the primordial one, with 
the debris of the destroyed clusters constituting a large 
fraction of the spheroidal (bulge + halo) stellar popula- 
tion. Mackey & van den Bergh (j2D05) by means of ob- 
servational differences in properties of three Galactic GC 
subsystems estimate that the present population could be 
a fraction ~ 2/3 of the original one. Mackey & Gilmore 
(2004) estimate a lower limit of 50% for the destruction 
rate over the last Hubble time, an intermediate value be- 
tween those of Gnedin & Ostriker (1992J) and Mackey & 
van den Bergh l|2()()5jl . 

At large Galactocentric distances the efficiency of the 
GC destruction mechanisms should be minimized with re- 
spect to the central region. This assumption is supported 
by numerical simulations showing that low-concentration, 
high-mass GCs are efficiently destroyed in the inner halo 
but are able to survive at Rgc > 10 kpc (Vesperini & 
Heggie 1997). For large Rgc the number-density of GCs 
in RS-GC05 is well described by the single power-law 
p(R) = (14.2 ± 3.3)R- (3 - 9±01 ), the discrepancy with re- 
spect to the observed profile becoming increasingly larger 
for Rqc < 3 kpc (Fig. 0- panel (c)). 

An estimate of the past destruction rate can be de- 
rived by comparing the number of GCs in RS-GC05 with 
that corresponding to the extrapolation of the large-Rcc 
power-law to the inner regions, under the assumption that 
the difference in both numbers is basically due to the cu- 
mulative past destruction of GCs. This estimate should 
be taken as lower-limit since we neglect GC destruction 
for large Rgc ■ To derive this value we integrate the large- 
Rcc power-law through the Galactocentric distance range 
over which the 116 GCs of RS-GC05 are observed. Taking 
into account the uncertainties in the parameters of the 
power-law fit we estimate that the present GC population 
represents a fraction of < 23 ± 6% of the primordial one. 
Thus, a lower-limit for the past destruction rate is ~ 77% 
over a Hubble time. This estimate is compatible with that 
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in Gnedin & Ostriker ( 1997) and about three times larger 
than that derived by Mackey & van den Bergh (|2005). 

The above estimate is based on the assumption that 
the flattening in the central radial profile is essentially 
due to GC destruction. However, if the flattening is partly 
primordial, as suggested by Parmentier & Grebel ( 2005 ) , 
the ~ 77% destruction rate should in fact be taken as an 
approximate upper limit. In this case our estimate would 
agree with the central range of Gnedin & Ostriker (1997). 

4.6. The central kpc 

With the present determination of the Galactocentric dis- 
tance (Sect. I3.1|) . 9 GCs of the reduced sample are lo- 
cated within 1 kpc of the Galactic center (H03 contains 11 
such GCs) . The metal-poor GCs arc Palomar 6 and HP 1 
at R G c < 0.5 kpc, and NGC6355, Terzan9, NGC6522, 
NGC6558, and NGC6401 at 0.5 < R G c(kpc) < 1.0. 
The metal-rich ones are Terzan5 and ES0456SC38 at 
0.6 < Rcc(kpc) < 0.7. Errors affect these individual esti- 
mates, but the ensemble might give hints on the extent of 
a potential avoidance zone, or a central region of enhanced 
destruction rates (e.g. Aguilar, Hut & Ostriker 1988). 

5. Discussion 

In the analysis of the reduced sample GCs we found that 
the radial density profiles of the metal-rich and metal- 
poor GCs are well described by a composed power-law of 
the form p(R) oc (1 + R/R c ) with 3.9 < a < 4.4. 
For Galactocentric distances larger than ~ 2 kpc both 
observed profiles fall off as a single power-law R~ n with 
3.2 < n < 3.9, while for inner regions they both flatten 
in a similar way. Structurally, the only difference besides 
geometry (Sect.OJ) between the spatial distributions of the 
metal-rich and metal-poor GCs is that the former extends 
basically to the Solar circle, while the latter spans from the 
central parts to the outer halo. This suggests that a signif- 
icant fraction of the metal-rich and metal-poor GCs share 
a common origin. These conclusions are not significantly 
affected by the oblate and nearly spherical geometries of 
the metal-rich and metal-poor sub-systems, respectively 
(Sect. El. 

As pointed out by Djorgovski & Meylan ( 1994), the an- 
alytical function that we adopted to fit the radial density 
profiles does not have a physically consistent counterpart. 
However, the steep slopes implied by this function for the 
GC radial density profiles rule out scenarios involving the 
pure monolithic collapse of an isothermal, uniform density 
cloud. This kind of collapse would produce flatter radial 
density profiles (Abadi, Navarro & Steinmetz 120051 and 
references therein). In addition, RS-GC05 radial density 
slopes are significantly steeper than that of the dark mat- 
ter halo (Merritt et al. l2005f) . which in principle precludes 
a common origin of these structures. Consequently, addi- 
tional mechanisms might have been necessary to increase 
the density of GCs in the central regions of the Milky Way, 
such as mergers in the early phases. 



The rather steep density distribution of the stellar 
halo was previously derived using globular clusters, RR 
Lyrae stars, blue horizontal-branch (BHB) stars, and star 
counts. Harris Ijl976(l and Zinn l|1985J) have shown that 
the metal-poor GCs distribute radially following a R _n 
power-law profile, with n = 3.5. Using observations of RR 
Lyrae Hawkins 119841 derived n = 3.1 with an axial ratio 
c/a = 0.9. Bahcall & Soneira lfT^54) and Gilmore iflWjl 
found c/a w 0.8. Using BHB stars Preston, Shectman & 
Beers (|1991|) found an increase in the axial ratio from 
c/a « 0.5 to c/a w 1 up to 20 kpc with n = 3.5. Recently, 
Yanny et al. (2000$ used BHB tracers from Sloan Digital 
Sky Survey data to derive c/a = 0.65 and n = 3.2. Ivezic 
et al. (|2000l) found that the RR Lyrae column density fol- 
lows a shallower power law with n = 2.7. Robin, Reyle & 
Crczc ( 2000 ) from deep wide field star counts estimated a 
halo flattening of 0.76, and n = 2.44. 

Evidence of merger events in galaxies has considerably 
increased in recent years, both on observational and the- 
oretical grounds. Deep observations of stars in luminous 
halos associated with numerical simulations of galaxy for- 
mation in A-cold dark matter (ACDM) scenarios indicate 
that a large fraction of the stellar content in the halo was 
not formed in situ. Abadi, Navarro & Steinmetz H2()()5l and 
references therein) suggest that this fraction may have 
been accreted from protogalaxies during earlier merger 
events. The resulting mass density profiles behave as a 
power-law R _n with n = 3 at the luminous edge of the 
galaxy and n > 4 more externally, at the virial radius. 
In addition, the density profile of the outer stellar halo is 
more centrally concentrated with a steeper slope than the 
dark matter halo, whose density profile is caracterized by 
slopes n = 1 — 3 (Merritt et al. 2005 ) . Hierarchical galaxy 
formation models under the ACDM framework that are 
successful in reproducing the radial density profile of the 
Milky Way stellar halo indicate that this structure formed 
from ~ 100 tidally disrupted, accreted dwarf satellites 
(Bullock, Kravtsov & WeinbergBOOI)) . What emerges from 
this is an observational/theoretical picture for the forma- 
tion and evolution of the stellar halo involving violent re- 
laxation and accretion, consistent with hierarchical mod- 
els of galaxy formation (Bellazzini, Ferraro & Ibata 2003 
Abadi, Navarro & Steinmetz I2005|l . 

In this sense, the observed similarity of the metal-rich 
and metal-poor radial density profiles of the Galactic GC 
system with that of the stellar halo suggests that, in ad- 
dition to the GCs formed in the primordial collapse, a 
non-negligible fraction of the GCs presently in the Milky 
Way was probably accreted during an early period of ac- 
tive merging. This scenario seems to apply to the bulge as 
well, which raises the possiblity of an early merger affect- 
ing the central parts of the Galaxy. In the present universe 
there are examples of such mergers, e.g. NGC 1275 (Zcpf 
et al. 113351 Holtzman et al. lT332"|) . 

Additional support for this scenario comes from the 
fact that the radial density profiles of the GCs of the 
reduced sample are equally well fitted by Sersic's law, 
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p(R) oc e -b[(R/Ro) (1/n) -i] ) with n K 41 ( mc tal-poor), 
n « 2.9 (metal-rich), and n « 3.0 (combined metal- 
poor/metal- rich GCs). Sersic's law with 2 < n < 4 is 
thought to apply to systems resulting from the mixing 
that follows from violent relaxation or merging (Merritt 
et al. 120051 and references therein). Besides, chemical evo- 
lution models that reproduce the observed abundances of 
stars in the bulge suggest that the Galactic bulge formed 
from the same gas but faster than the inner Galactic 
halo (Matteucci & Romano 119991 Matteucci, Romano & 
Molaro 1999). The minimum at about [Fe/H] w -0.75 in 
the observed metallicity distribution of GCs (Fig. may 
reflect an external mechanism such as merging to explain 
the exceedingly large number of metal-rich GCs. An early 
merger in the Milky Way with a relatively massive galaxy 
might have provided the excess metal-rich star formation 
in the central parts. 

Further evidence on the bulge formation via collapse 
and/or additional mechanisms will be given by detailed 
derivation of metallicities and abundance ratios in com- 
prehensive samples of GCs. However, such detailed infor- 
mation for bulge GCs is presently scarce. 

6. Concluding remarks 

In this paper the Galactic globular cluster system was de- 
contaminated of the objects with strong evidence of exter- 
nal origin and/or ages younger than the Galaxy collapse. 
The resulting reduced sample contains 116 GCs, 81 metal- 
poor and intermediate metallicity clusters ([Fe/H] < 
—0.75), 33 metal-rich and 2 with unknown metallicity. 
The classical bimodal metallicity distribution is enhanced 
in the reduced sample. 

Projections of the observed heliocentric distances onto 
the (x,y,z) planes show that the metal-rich GCs distribute 
in a central region of dimensions ~ 12 kpc x 1 1 kpc x 5 kpc, 
whose structure resembles an oblate spheroidal with axial 
ratio c/a w 0.4. The metal-poor ones span a region of di- 
mensions ~ 35 kpc x 36 kpc x 30 kpc, with a shape similar 
to a slightly flattened sphere with c/a ~ 0.8. The metal- 
poor GCs in the reduced sample extend into the beginning 
of the outer halo. Based on the projected number-density 
of GCs along the x-direction we measured the distance 
of the Sun to the Galactic center as Ro = 7.2 ± 0.3 kpc. 
This value was obtained considering the spatial distribu- 
tion of 116 GCs, and is ~ 10% smaller than the widely 
used estimate of Reid (|1993[l . 

Based on structural similarities of the radial density 
profiles of the present-day GC population with the stel- 
lar halo one can build a scenario where, besides the GCs 
formed in the primordial collapse, a non-negligible frac- 
tion of the Milky Way GCs was probably accreted from 
satellites during an early period of merging. Observational 
and theoretical evidence support this picture, e.g. the GCs 
formed as a consequence of mergers in NGC 1275 (Zepf et 
al. lT§9"5l Holtzman et al. IT9"9"2"|l . 

The present decontamination procedure was not sen- 
sitive to all accretions that may have occurred in the first 



Gyrs of the Galaxy, including an eventual major merging, 
since the observed radial density profiles still appear to 
preserve traces of the earliest merger (s). 

Assuming that the flattening in the observed radial 
density profiles is a consequence of the cumulative GC de- 
pletion mostly by bulge and disk shocking we estimated 
that the present GC population represents a fraction of 
< 23 ±6% of the original one. This in turn implies a lower- 
limit destruction rate of ~ 77% over a Hubble time. This 
estimate is compatible with that of Gnedin & Ostriker 
( 1997| and somewhat larger than that derived by Mackey 
& van den Bergh (2005). However, if the central flatten- 
ing is partly primordial (Parmentier & Grebel 2005) our 
estimate would in fact be an upper limit. 

The significant improvement in the accuracy of GC 
data over the last years, as analysed in the present work, 
has shed light on the issue whether the metal-rich GCs are 
associated to a disk (e.g. Zinn 1 198*51 Armandroff 1985J) or 
a spheroidal subsystem. The fact that the volume-density 
radial distribution of GCs of the reduced sample can be 
described both by a core-like power-law or a Sersic's law 
indicates that the metal-rich GC subsytem is spheroidal. 

The present study pointed out that besides the GC 
accretions from dwarfs and/or formation later than the 
primordial collapse, the radial density distributions re- 
quire, in addition to a primordial collapse component 
a non-negligible early merger population. This scenario 
provides also a natural explanation to the second peak 
in the bimodal metallicity distribution. Through gravita- 
tional lensing, large galaxies at high-redshift have been 
detected in the starburst stage, in an epoch compatible 
to that of the Milky Way's primordial collapse. Examples 
are the z ~ 7 and M ~ 10 9 M© galaxy lensed by the 
Abell2218 cluster (Egami et al. I2005JI . and the z ~ 5.5, 
M ~ 1 — 6 x 1O 1O M0 starforming galaxy in the field of the 
cluster RDCS 1252.9-2927 (Dow-Hygelund et al. EES). 
Collapse or its combination with merging are supported by 
the Galactic GCs and large redshift observations of galax- 
ies. On the other hand pure hierarchical galaxy formation 
has yet to be observed in detail. 
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Table 1. Updated Globular Cluster parameters - the GC05 sample 
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.43 


30.7 


30.0 


30.0 


6.59 


22.46 


19.86 




GC1-47 


Pall3 


87 


10 


-42.70 


0.05 


-1 


.74 


26.9 


27.8 


27.6 


1.00 


19.74 


-18.24 


R 7 


Pegasus, GC1-124 


Pall 


130 


07 


19.03 


0.15 


-0 


.60 


10.9 


17.0 


16.3 


-6.63 


7.89 


3.55 


Y 3 CT 9 


GC1-6 


NGC288 


152 


28 


-89.38 


0.03 


-1 


24 


8.3 


11.6 


11.0 


-0.08 


0.04 


-8.30 


R 8 


GC1-2.ES0474SC37 


SDSSJ1049+5103+ 


155. 


65 


55.62 


0.01** 


-1.70* 


45.0 


50.0 


49.1 


-23.15 


10.48 


37.14 






Whiting 1* 


161 


62 


-60.64 


0.04 


-1.20 


45.0 




48.1 


-20.94 


6.96 


-39.22 


Y 4 


WHI,B0200-03 


Pal2 


170 


53 


-9.07 


1.24 


-1.30 


27.6 


35.4 


34.7 


-26.88 


4.48 


-4.35 




GC1-7 


NGC2419 


180 


37 


25.24 


0.11 


-2.12 


84.2 


91.5 


90.8 


-76.16 


-0.49 


35.90 


DN 2 


GC1-12 
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Table 1. Continued 



E. Bica et al.: Galactic system of globular clusters 



Heliocentric 



GC 


t 


b 


E(B-V) 


[Fc/H] 


d 


d GC 


d GC 


X 


y 


z 


Notes 


Alternative designations 




(°) 


(°) 










(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Pal4 


202.31 


71.80 





01 


-1 


.48 


109.2 


111 


8 


111 


5 


-31.55 


-12.95 


103.74 




Ursa Majoris,GCl-17 


Eridanus 


218.11 


-41.33 





.02 


-1 


.40 


90.2 


95. 


2 


94 


.6 


-53.29 


-41.80 


-59.57 


yll 


ESO551SCl,C0422-213 


M79 


227.23 


-29.35 


(1 


01 


-1 


.57 


12.9 


18. 


8 


18. 


.1 


-7.64 


-8.25 


-6.32 


C 9 


NGC1904,GCl-10,ESO487SC7 


Pal3 


240.14 


41.86 





04 


-1 


.66 


92.7 


95. 


9 


95 


6 


-34.37 


-59.88 


61.86 


yll 


SextansC 


NGC1851 


244.51 


-35.04 





(12 


-1 


.22 


12.1 


10. 


7 


16 


.1 


-4.26 


-8.94 


-6.95 




GCl-9,ESO305SC16 


NGC2298 


245.63 


-16.01 





14 


-1 


.85 


10.7 


15 


7 


15 


,1 


-4.24 


-9.37 


-2.95 


c 9 


GC1-11,ES0366SC22 


NGC4147 


252.85 


77.19 





02 


-1 


.83 


19.3 


21 


3 


21. 





-1.26 


-4.09 


18.82 


R 5 ST 9 


GC1-18 


AMI 


258.36 


-48.47 





00 


-1 


.80 


121.9 


123 


2 


123 


1 


-16.31 


-79.16 


-91.26 




El,ESO201SC10 


Pyxis 


261.32 


7.00 





21 


-1 


.30 


39.4 


41. 


.4 


41 


1 


-5.90 


-38.66 


4.80 




Wcinbcrgcr3 


NGC1261 


270.54 


-52.13 





.01 


-1 


.35 


16.4 


18. 


2 


17. 


.9 


0.09 


-10.07 


-12.95 


yll 


GC1-5,ES0155SC11 


NGC3201 


277.23 


8.64 





.23 


-1 


.58 


5.2 


9 





8 


,4 


0.65 


-5.10 


0.78 




GC1-15,ES0263SC26 


NGC2808 


282.19 


-11.25 


o. 


22 


_1 


.15 


9.3 


11 





10 


.5 


1.93 


-8.92 


-1.81 


C 9 Y 11 


GC1-13,ES091SC1 


E3 


292.27 


-19.02 


o 


.30 


-0 


.80 


4.3 


7 


(j 


6 


.9 


1.54 


-3.76 


-1.40 




ES037SC1 


M68 


299.63 


36.05 


o. 


.05 


-2 


.00 


10.2 


10 


1 


9 


9 


4.08 


-7.17 


6.00 




NGC4590,GCl-20,ESO506SC30 


Rul06 


300.89 


11.67 


o. 


2(1 


_1 


.07 


21.2 


18. 


5 


18 


g 


10.66 


-17.82 


4.29 


Y 3 


OCl-887,ESO218SC10 


NGC4372 


300.99 


-9.88 


() 


.39 


_2 


.09 


5.8 


7 


1 


6 


.6 


2.94 


-4.90 


-1.00 




GC1-19,ES064SC6 


NGC362 


301.53 


-46.25 


Jj 


''9 






8.5 








! ' 


3.07 


-5.01 


-6.14 


Y^R 12 


GC1-3 ES051SC13 


NGC4833 


303.61 


-8.01 




' 




t r 


0.0 




' 






3.29 


-4.95 


-0.84 




GC1-21,ES065SC4 


47 Tucanac 


305.90 


-44.89 










4.5 










1.87 


-2.58 


-3.18 




NGC104,GCl-l,ESO50SC9 


IC4499 


307.35 


-20.47 


0. 


.23 


- 1 


.00 


18.9 


1 5 . 


7 


15 


11 


10.74 


-14.08 


-6.61 


Y 3 


GC1-30 ES022SC5 


u Centauri 


309. 10 


14.97 





12 


-1 


.02 


5 . 3 


6 


.4 


5 


.8 


3.23 


-3.97 


1.37 




IN (j(^o toy ,001-z4,tjDOi (UbOll 


NGC5286 


311.61 


10.57 





.24 


-1 


.07 


11.0 


8 


.4 


8 


.3 


7.18 


-8.08 


2.02 


CT 9 


GCl-26,ESO220SC38 


NGC6101 


317.75 


-15.82 





05 


-1 


.82 


15.3 


11 


1 


11 


.3 


10.90 


-9.90 


-4.17 




GCl-40,ESO69SC4 


AM4 


320.28 


33.51 





04 


-2 


.00 


29.9 


25 


5 


25 


.9 


19.18 


-15.93 


16.51 




AM1353-265 


NGC6362 


325.55 


-17.57 





0!) 


-0 


.95 


8.1 


5 


.3 


5 


1 


6.37 


-4.37 


-2.45 




GCl-66,ESO102SC8 


NGC5927 


326.60 


4.86 





45 


-0 


.37 


7.6 


4 


.5 


4. 


3 


6.32 


-4.17 


0.64 




GC1-35,ES0224SC4 


IN LrO094D 


327.58 


4. 19 





54 


-1 


.38 


12.8 


7 


4 


7 


.8 


10.78 


-6.84 


0.94 




104DD0,GO-3b,Eb(J224b0 / 


BH176 


328.41 


4.34 





.77 


-0 


.13 


14.5 


8 


.8 


9. 


.2 


12.32 


-7.57 


1.10 


^10 


ES0224SC8 


Lynga7 


328.77 


-2.79 





.73 


-0 


.62 


7.2 


4 


2 


3 


.9 


6.15 


-3.73 


-0.35 




OCl-949,BH184,ESO!78SCll 


NGC5694 


331.06 


30.36 





(1!) 


-1 


.86 


34.7 


29 


1 


29 


(i 


26.20 


-14.49 


17.54 




GCl-29,ESO512SC10 


NGC5824 


332.55 


22.07 





.13 


-1 


.85 


32.0 


25. 


8 


20 


4 


26.32 


-13.67 


12.02 




GC1-31,ES0387SC1 


M53 


332.96 


79.76 





02 


-1 


.99 


18.3 


18. 


8 


18 


(i 


2.90 


-1.48 


18.01 




NGC5024,GCl-22 


NGC5053 


335.69 


78.94 





04 


-2 


.29 


16.4 


10. 


9 


Hi 


.7 


2.87 


-1.30 


16.10 




0267,00-970, GC1-23 


NGC6752 


336.50 


-25.63 





04 


-1 


.56 


4.0 


5 


.2 


4 


.5 


3.31 


-1.44 


-1.73 




GO-108,ESO141SC30 


NGC5986 


337.02 


13.27 





28 


-1 


.58 


10.5 


4 


.8 


5 


.1 


9.41 


-3.99 


2.41 




G0-37,ES0329SC18 


NGC6397 


338.17 


-11.96 





.18 


-1 


.95 


2.3 


(i 





5. 


.2 


2.09 


-0.84 


-0.48 




G0-74,ES0181SC4 


NGC6352 


341.42 


-7.17 





.21 


-0 


.70 


5.7 


3 


3 


2. 


.7 


5.36 


-1.80 


-0.71 


Y 11 


Mel-170,Cr328, 00-993, G0-64,ES0228SC3 


NGC6584 


342.14 


-16.41 





10 


-1 


.49 


13.4 


7. 





7. 


4 


12.23 


-3.94 


-3.79 




G01-92,ESO229SC14 


NGC5634 


342.21 


49.26 





05 


-1 


.88 


25.9 


21 


9 


22 


.1 


16.09 


-5.16 


19.62 




GC1-28 


NGC6139 


342.37 


6.94 





75 


-1 


.68 


10.1 


3 


(i 


4. 





9.56 


-3.04 


1.22 




G01-43,ESO331SC4 


NGC5897 


342.95 


30.29 





(111 


-1 


.80 


12.8 


7 


.7 


8 





10.57 


-3.24 


6.46 




GC1-33,ES0582SC2 


Tz3 


345.08 


9.19 





72 


-0 


.73 


7.5 


2 


.4 


2. 


.2 


7.15 


-1.91 


1.20 




ESO390SC6 


NGC6388 


345.56 


-6.74 





.37 


-0 


.00 


11.5 


4 


4 


5 





11.06 


-2.85 


-1.35 




GCl-70,ESO279SC2 


ESO280SC6 


346.90 


-12.57 





07 


-2 


.00 


21.7 


14 


5 


15. 





20.63 


-4.80 


-4.72 






NGC6256 


347.79 


3.31 


1 


03 


-0 


.70 


6.6 


2 


1 


1 


6 


6.44 


-1.39 


0.38 




BH208,ESO391SC6 


NGC6496 


348.02 


-10.01 





.15 


-0 


.04 


11.5 


4 


.3 


4. 


9 


11.08 


-2.35 


-2.00 




GCl-80,ESO279SC13 


NGC6541 


349.29 


-11.18 





14 


-1 


.83 


7.0 


2 


2 


1. 


.9 


6.75 


-1.28 


-1.36 




GCl-86,ESO280SC4 


NGC6380 


350.18 


-3.42 


1. 


.17 


-0 


.50 


10.7 


3 


.2 


3 


.8 


10.52 


-1.82 


-0.64 




Tonl,Pis25,GO-68,BH233,ES0333SC14 


Ton2 


350.80 


-3.42 


1 


24 


-0 


.50 


8.1 


1 


.4 


1. 


(i 


7.98 


-1.29 


-0.48 




Pis26,GO-71,BH236,ES0333SC16 


M4 


350.97 


15.97 





36 


-1 


.20 


2.2 


5 


.9 


5 


.2 


2.09 


-0.33 


0.61 




NGC6121,GC1-41,ES0517SC1 


ES0452SC11 


351.91 


12.1 





40 


-1 


.50 


7.8 


2 





2 





7.55 


-1.07 


1.64 




C1636-283 


NGC6144 


351.93 


15.70 





36 


-1 


.75 


10.3 


3 


6 


4 


.1 


9.82 


-1.39 


2.79 




G0-42,ES0517SC6 


M80 


352.67 


19.46 





.18 


-1 


.75 


10.0 


3 


.8 


4 


.1 


9.35 


-1.20 


3.33 




NGC6093,GCl-39,ESO516SCll 


NGC6441 


353.53 


-5.01 





47 


-0 


.53 


11.2 


3 


5 


4 


.2 


11.09 


-1.26 


-0.98 




GC1-78,ES0393SC34 


M62 


353.58 


7.32 





47 


-1 


.29 


6.9 


1 


.7 


1. 


.2 


6.80 


-0.77 


0.88 




NGC6266,GC1-51,ES0453SC14 


Lillcrl 


354.84 


-0.16 


3 


06 


0.22 


10.5 


2 


(i 


3 


4 


10.46 


-0.94 


-0.03 






NGC6453 


355.72 


-3.87 





66 


-1.53 


11.2 


3 


.3 


4 


.1 


11.14 


-0.83 


-0.76 




GC1-79,ES0393SC36 


NGC6304 


355.83 


5.38 





.53 


-0.59 


0.1 


2 


.1 


1. 


.4 


6.06 


-0.44 


0.57 




GC1-56,ES0454SC2 


Tz4 


356.02 


1.31 


2 


.35 


-1.60 


9.1 


1. 


.3 


2 





9.08 


-0.63 


0.21 




HP4,ES0454SC7 


Tz2 


356.32 


2.30 


1 


57 


-0.40 


8.7 





.9 


1 


6 


8.68 


-0.56 


0.35 




HP3,BII228.ES()454SC29 


Djorgovskil 


356.67 


-2.48 


1. 


44 


-2.00 


8.8 


1 





1. 


.7 


8.78 


-0.51 


-0.38 






M19 


356.87 


9.38 





41 


-1.68 


8.7 


1 


6 


2 





8.57 


-0.47 


1.42 




NGC6273,GC1-52,ES0518SC7 


NGC6316 


357.18 


5.76 





.51 


-0.55 


11.0 


3 


2 


3 


.9 


10.93 


-0.54 


1.10 




G0-57,ES0454SC4 


HP1 


357.42 


2.12 





74 


-1.55 


7.4 





.8 





.5 


7.39 


-0.33 


0.27 




GC1-67,BH229,ES0455SC11 


Tzl 


357.57 


1.00 


2 


28 


-1.30 


6.2 


1 


.8 


1 


.1 


6.19 


-0.26 


0.11 




HP2,GC1-69,ES0455SC23 


NGC6293 


357.62 


7.83 





41 


-1.92 


8.8 


1 


4 


1. 


.9 


8.71 


-0.36 


1.20 




GC1-55,ES0519SC5 


NGC6284 


358.35 


9.94 





28 


-1.32 


14.7 


6 


.9 


7 


7 


14.47 


-0.42 


2.54 




G0-53,ES0518SC9 


Tz6 


358.57 


-2.16 


2 


14 


-0.50 


9.5 


1 


.6 


2 


.3 


9.49 


-0.24 


-0.36 




HP5,BH249,ESO393SC36,BDSB103 


NGC6235 


358.92 


13.52 





.36 


-1.40 


10.0 


2 


.9 


3 


.4 


9.72 


-0.18 


2.34 




GC1-58,ES0519SC11 


NGC6355 


359.58 


5.43 





75 


-1.50 


7.2 


1 








7 


7.17 


-0.05 


0.68 




O330, 00-1036, G0-63,ES0519SC15 



Table Notes. S: Sagittarius; C: Canis Major; Y: Younger than lOGyr; R: Retrograde orbit; DN: Dwarf nucleus; ST: Sagittarius 
tidal tail; CT: Canis Major tidal tail. (1) - da Costa & Armandroff H1995H : (2) - Mackey & van den Bergh (120051 : (3) - 
Rosenberg et al. ftMty : (4) - Carraro SftWty : (5) - van den B ergh ftiEEty : (6) - Cote e t al. fT335) : (7) - Siegel et al. pi) ; 
(8) - Guo et al. l(T993jl : (9) - Forbes, S trader & Brodie p)04)l : (10) - Phelps & Schick pOOSjl : (11) - Salaris & Weiss <|2002jl ; 
(12) - Dinescu, Girard & van Altena (1999). (f) and (\) cluster parameters from Willman, Blanton & West et al. (2005) 
and Carraro I2005H . respectively. (*): based on the conclusion by Willman, Blanton & West et al. (2005) that the cluster is 
more metal-poor than Pal 5; (**): based on Schlegel, Finkbeiner & Davis 1)19980 . Col. 7: Galactocentric distance from H03. 
Col. 8 : this work. 



